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Volcanic SO2 clouds measured by OUS 
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TOMS: Total Ozone 
Mapping Spectrometer 

1 978-2005 


Low sensitivity 
Low spatial resolution 
No altitude information 
No ‘monitoring’ capability 



Pre-Aura view of global volcanic SO2 emissions 
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A time-averaged inventory of subaerial volcanic 
sulfur emissions 

RJ, Andres and A,D. ICasgnoc 

InsliUiti? of Ncnthtm Engineering, University of Alaska Fairl>anks 


Abstract. A time-averaged inventory of subaerial volcanic sulfur (S) emissions was 
compiled primarily for the use of global S and sulfate modelers. This inventory relies 
upon the 25-year history of S, primarily sulfur dioxide (SO 2 ), measurements at 
volcanoes. Subaerial volcanic SO^ emissions indicate a 13 Tg/a SOj time-averaged 
flux, based upon an early 1970s to 1997 time frame. When considering other S species 
present in volcanic emissions, a time-averaged inventory of subaeriai volcanic S fluxes 
is 10,4 Tg/a S. These time-averaged fluxes are conservative minimum fluxes since 
they rely upon actual measurements. The temporal, spatial, and chemical 
inhomogeneities inherent to this system gave higher S fluxes in specific years. Despite 
its relatively small proportion in the atmospheric S cycle, the temporal and spatial 
distribution of volcanic S emissions provide disproportionate effects at local, regional, 
and global scales. This work contribute.? to the Global Emissions Inventory Activity. 
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• Volcanic S emission inventories ‘static’ and 
used ground-based data 


Anthropogenic (current) 



Motivation: 

• Climate impact of volcanic emissions 

• Global fluxes of other volatile species 

• Field sites for volcanic gas studies 



[Grafetal., 1997; Andres & Kasgnoc, 1998; Halmeretal., 2002; Smith etal., 2011] 








OMI monitors SOo emissions in South^America 
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Average SO 2 map from OMI 
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Cam, S.A., N.A. Krotkov, A.J. 
Kmeger, K. Yang, RF. Levelt, 

"Sulfur dioxide emissions from 
Pemvian copper smelters detected by 
the OMI", GRL, 2007. 


• SO 2 emissions in Ecuador and Colombia from Tungurahua, 
Reventador, Galeras and Huila volcanoes 

• Volcanic emissions from Ubinas volcano (Peru) drift over No region 
in 2006-2007 

• Upgraded sulfur capture technology at No reduces SO 2 emissions 

• Shut-down of La Oroya smelter in mid-2009; reduction expected 





Volcanic SO2 emissions inventories are inaccurate 
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• OMI measurements indicate 
deficiencies in current ‘bottom 
up’ volcanic SOg emission 
inventories used in CTMs 


REMOTE model simulation of annual 
mean SO2 columns over Indonesia 
[Pfefferetal.,ACP, 2006] 
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SOg absorption coefficient [cm'^] 


Identification of global volcanic SO2 sources 
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• SO 2 index defined using OMI ozone algorithm residuals (sensitive to SO 2 ) 

• Rank top 1000 values of SO 2 index on each day (lat 65-S - 80-N) 

• Locate active volcanoes within 50 km radius of these OMI pixels (using modified 
Smithsonian Global Volcanism Program database) 

• Record unique volcanic sources for each day 

• Repeat for each day of OMI measurements (>3200 from Sep 2004 - Dec 2013) 





Rank 


Persistent volcanic SOg sources (2004-2013; ~3200 days) 
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Infrequently monitored 
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[Andres and Kasgnoc, 1998] _ 
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Most persistent volcanic SO 2 sources in 2013 
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SO2 emission rate estimation from satellite data 
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M = SO2 mass in pixel (kg) 

V = wind speed (m s '*) 

L = length of plume (m) 

Q = SO2 flux (kg s'^) 

• Similar approach used to estimate 
smoke and NO 2 emissions from fires 
[Ichoku and Kaufman, 2005; Mebust 
et al., 2011] 

• Note that asymmetry of OMI pixel 
affects plume detection 

• Chemistry correction [Mebust et al., 
2011] can be applied if SO 2 lifetime is 
known 


[Cam etal., Spec. Pub. Geol. Soc. Lon., 380, 2013] 
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• New volcanic SO2 emissions - 

database for CTMs - 

• Jun Wang’s talk (Tuesday) ~ 
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[Oppenheimer et a!., 2011] 
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Validation with ground-based SO2 measurements 



Ground-based SO2 flux (tons/day) 

• Good agreement with long-term, ground-based datasets 




NASA^KNMI/NI VR/FMI Contact Simon Cam (scarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/20/2006 17:00-18:41 UT 

Mass; 135.1 33 kt; Area; 202457 km^ SO^ max; 146.85 DU a1 Ion; -64.79 lal; 1 5.72 


'^0.2 Tg SO 2 



Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 




NASA''KNMI/NI VR'^FMI Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/21/2006 17:40-19:25 UT 

Mass: 195.725 k1; Area: 1400753 SO^ max: 32.71 DU al Ion: -74.40 lal: 1 1.42 



Normalised SOg column 
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[Cam et al., 2007; Praia et al., 2007] 



NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/22/2006 16:48-20:08 UT 

Mass: 188.536 k1; Area: 1880299 SO^ max: 23.66 DU al Ion: ^1.33la1: 14.47 
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[Cam et al., 2007; Praia et al., 2007 ] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 



Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 
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Mass: 171.187 k1; Area: 2410030 km^ SO^ max: 19.89 DU al Ion: -93.38 1 al: 13.66 


Soufriere Hills volcano (Montserrat) eruption, May 2006 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/24/2006 16:33-21 :35 UT 

Mass: 159.341 k1; Area: 2532097 SO^ max: 20. 18 DU alien: -107.02 lal: 12.14 



Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/25/2006 14:03-22:18 UT 

Mass: 168.052 k1; Area: 3126612 SO^ max: 18.41 DU al Ion: -114.18 lal: 13.47 



[Cam et al., 2007; Praia et al., 2007 ] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/26/2006 16:19-23:02 UT 

Mass: 146.936 k1; Area: 3468046 SO^ max: 12.14 DU alien: -121.13 lal: 13.78 



Normalised SOg column 



[Cam et a!., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/27/2006 20:21-23:45 UT 

Mass: 110.209 k1; Area: 2575001 SO^ max: 8.59 DU al Ion: -128.96 1 a1: 14.21 



Normalised SOg column 
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[Cam et al., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/28/2006 00:22-22:50 UT 

Mass: 108.896 k1; Area: 2279174 SO^ max: 8.1 1 DU al Ion: -135.54 1 a1: 13.82 



Normalised SOg column 
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[Cam et al., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/29/2006 01 :07-23:33 UT 

Mass : 98.493 kt; Area: 1701411 SO^ max: 6.22 DU a1 Ion: -145.52 lal: 12.46 



Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/30/2006 00:11-22:37 UT 

Mass: 107.076 k1; Area: 2034725 SO^ max: 6.74 DU al Ion: -146.34 lal: 1 3.83 



Normalised SOg column 
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[Cam et al., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 05/31/2006 00:54-23:21 UT 

M ass : 96.573 kt ; Are a: 2063998 km^SO^ max: 6. 10 DU a1 Ion: -156.59 lal: 11.76 



Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 06/02/2006 00:00-24:00 UT 

Mass: 92.692 kt; Area: 2247222 SO^ max: 5.54 DU a1 Ion: -165.01 lal: 9.94 



Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 06/03/2006 00:41-04:04 UT 

Mass: 82.835 kt; Area: 22D9377 SO^ max: 3.85 DU a1 Ion: -175.63 lal: 9.17 




Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 

40 50 60 

70 

80 

90 

100 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 06/04/2006 01 :25-06:23 UT 

Mass: 79.103 kt; Are a: 2352265 SOj max: 3.01 DU a1 Ion: 155.20 lal: 1 1.69 




Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 
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NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 


Soufriere Hills volcano (Montserrat) eruption, May 2006 


Aura/OMI - 06/05/2006 02:08-07:05 UT 

Mass: 71.854 kt; Are a: 2 170677 SOj max: 2.94 DU a1 Ion: 151.98 lal:8.57 




Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 
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NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 



Normalised SOg column 
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[Cam et a!., 2007; Praia et al., 2007] 
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Mass: 54.517 kt; Area: 1986410 SOj max: 2.67 DU al Ion: 138.03 lal:5.65 


Soufriere Hills volcano (Montserrat) eruption, May 2006 




NASA-''KNMI/NI VFl''FIVII Contact: Simon Cam (sGarn@umbc.edu) 
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Mass:45.728kt;Area: 1734106km^SOjmax:2.12DU a1 Ion: 126.1 3 lal: 4.70 


Normalised SOg column 
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No volcanic SO 2 cloud of similar size tracked for more than 7 days by TOMS 


Soufriere Hills volcano (Montserrat) eruption, May 2006 




New insights into volcanic SO2 plume.dispersion 
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• Improve sulfur chemistry I 
scheme in CTMs 
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Days after July 31, 2008 


[Krotkov et al., JGR, 2010; 
Wang et al., AGP, 201 3] 
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A-Train data for February 2014 Kelut eruption 





Color Enhanced Infrared Imagery (11p,m) 


Aqua MODIS 02/13/2014 - 18:10 UTG 
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MODIS image courtesy of NOAA/CIMSS 




A-Train data for February 2014 Kelut eruption 
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Kristiansen et al., GRL, in prep. 





Direct retrieval of SO2 altitude from UV radiances 
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• SO 2 altitude directly retrieved 
from hyperspectral UV radiances 

• Validate with CALIPSO, MLS 

• Critical for climate impact, 
aviation hazards and CTM source 
term 


[Yang et a!., JGR, 2010] 


SO 2 altitude retrievals for 2008 Kasatochi eruption 




UV satellite volcanic SO 2 inventory (1978 - 2014) 


333 eruptions; 94 Tg total SO 2 ; mean 0.28 Tg; 1 a = 1 .3 Tg 
~6 eruptions/yr (1 979-2004); ~1 9 eruptions/yr (2004-201 4) “ 
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■ TOMS data gap 
• Total explosive SO2 emission 
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Stratospheric AOD and the global warming ‘hiatus’ 
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Vernier et al., GRL, 2011 
Solomon et al., Science, 2011 
Santeretal., Nature Geosci., 2014 
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Increase in tropical stratospheric AOD since 2000 
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• Eruptions filtered using 
Volcanic Explosivity Index 
(VEI) 

• VEI not a good proxy for 
climate impact (i.e., SO2 
emissions) 

• VEI no longer needed -> 
use satellite observations 


Ulawun(UI> 

Shiveluch CSh) 

Ruang CRu) 

Reventador (Ra) 

Anatahan (At) 

Manam (Ma) 

Sierra Negra (Si) 

Soufriere Hills (So) 


Tavurvur (Ta) 

Jebel at Jair (Jb) 
Chaiten (Ch) 
Okmok (Ok) 
Kasatochi (Ka) 
Sarychev CSa) 
Eyjafjallajokull (Ey) 


Merapi (Me) 

Nabro (Na) 

Observed SAOD 

CERES net dear-sky SW (raw) 

CERES net dear-sky SW (filtered) 

RSSTLT(ENSOout) 

RSSTLT (ENSO out, filtered) 


Santer et al., Nature Geosci. [2014] 




SO 2 mass [kt] 


Volcanic SO 2 and stratospheric aerosol since 2002 
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Cam et al., in prep. 
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Stratospheric Aerosol Optical Depth at 550 nm 




OMI SO2 validation - Okmok (Aleutian Is) eruption 


Aura/OMI - 07/17/2008 00:00-24:00 UT 

SO 2 mass: 67.350 kl; Area: 968738 km^ SOj max: 23.27 DU al Ion: -143.54 lal: 40.71 ; 22:1 9UTC 
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MFDOAS-OMI comparison 
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Extreme SO2 columns in volcanic clouds 
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• SO 2 in fresh eruption clouds - highest trace 
gas columns measured 

• Challenge to validate extreme SO 2 column 
amounts - UAVs, balloons? 

• More vole. SO 2 validation needed in general 

• Complex radiative transfer (ash, hydrometeors) 



[Yang et al., GRL, 2009] 
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Sierra Negra - October 2005 
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Current Bardabunga/Holuhraun eruption (Iceland) 


Aura/OMI -09/05/2014 13:11-13:16 UT - Orbit 53943 

SO 2 mass: 19.60 kt; Area: 287562 km^; SO 2 max: 35.04 DU at Ion: -13.15 lat: 59.12 ; 13:13UTC 


Aura/OMI -09/16/2014 12:55-12:57 UT - Orbit 54103 

SO 2 mass: 0.24 kt; Area: 13502 km^; SO 2 max: 2.05 DU at Ion: -14.81 lat: 65.57 ; 12:56UTC 



SO 2 column TRM [DU] 



0.0 0.8 1.6 2.4 3.2 4.0 4.8 5.6 6.4 7.2 8.0 


SO 2 column TRM [DU] 



0.50 0.65 0.80 0.95 1.10 1.25 1.40 1.55 1.70 1.85 2.00 


Request from Icelandic Meteorological Office: We are getting very nice OMI images of the 
SO2 rich eruption cloud from the Holuhraun eruption. Could you help us to set up a way of 
turning the OMI observations into daily SO 2 mass and daily SO 2 emission rates? It's possible 
this eruption will last weeks or months and so we need such a technique to be automated 
and fast. ’ 


OMI SO2 data used in: U.S., Mexico, Guatemala, El Salvador, Costa Rica, Montserrat, Ecuador, Colombia, Peru, 
Chile, New Zealand, Vanuatu, Indonesia, Papua New Guinea, Philippines, DR Congo, Ethiopia, Russia, France 




Summary 


• Science contributions and applications of Aura/OMI SO2 

- Identification of new volcanic and anthropogenic SO2 sources 

- Updated voicanic SO2 emissions inventory 

- Effect of smali voicanic eruptions on ciimate 

- Radiative forcing of voicanic suifate aerosoi 

- Voicanic aerosol enhancement of lightning activity 

- Cioud seeding by volcanic aerosoi downwind of voicanoes 

- Lifetime of SO2 in upper troposphere and iower stratosphere 

- Suifur gas scavenging in eruption columns 

- Improving suifur chemistry in CTMs 

- Voicanic plume tracking for aviation hazard mitigation 

- Detection of eruptions in remote regions 

• Future developments in satellite SO2 observations 

- Improved spatial resolution (e.g., TROPOMI) -> better monitoring 
capability 

- Direct assimilation of volcanic SO2 observations into CTMs 




Thanks to: 

the Aura Science and OMI SIPS Teams 

NASA for support from the Aura Validation program, Aura Science Team, 

ACMAP program, MEaSUREs program 


Ambrym volcano (Vanuatu) - one of the strongest SO2 sources on Earth 
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